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Chapter 1 

Background of the Invention 



1.1 Field of the Invention 

Methods, machines and articles of manufacture for simulating financial vari- 
ables, calculating quantities, preparing to and/or executing transactions. 
US Classification 705-35 and/or related categories. 

1.2 Description of the Related Art 

Note that typing in the names of the authors or articles below in a search 
engine will in many cases allow access to the articles cited herein or related 
ones or ones by others. JSTOR allows access to many of the historical pa- 
pers. Individual subscriptions are available to many of the journals with 
access to JSTOR for that journal. Typing in several authors names together 

19 
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will pick up at least one site that maintains a bibliography of finance pa- 
pers in these areas. The Society of Actuaries Library, American Academy 
of Actuaries, Canadian Institute of Actuaries and Faculty and Institute of 
Actuaries contain works of the author and Mathematical Finance Company 
and many of these are available on-line at their websites. See the bibliogra- 
phy of term structure literature and derivative pricing by Don Chance [:>i] 
and the references in Duffie [15] and Duffie and Singleton [48], 

2-1 Pure Regime Switching 

See Bharucha-Reid's [21] Chapter 2 for a discussion of discrete space pro- 
cesses in continuous time, and Chapter 1 for processes discrete in time and 
space. See also Shiryaev [110] for regime switching in discrete time, i.e. 
Markov chains. The pure regime switching continuous time process implies 
a pure regime switching discrete time process, even with unequal time inter- 
vals. Moreover, matrix methods can be used to compute this matrix. See the 
references on matrix methods and other references cited elsewhere. See page 
64 of Bharucha-Reid [21]. For a finite time interval in which the generator 
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matrix is A, the transition matrix is given by 

oo 

P(t) = ^A n t n /n\ (1.1) 

This can be evaluated using the matrix methods discussed in the references. 
See particularly Patel, Laub and Doren [106], 

2-2 Regime Switching in Finance 

See Hamilton [GO] [CI] and [62] for use of regime switching in economics and 
finance. See Hamilton and Raj [03] for recent applications. Krolzig. See Ang 
and Bekaert [5]. Babbs and Nowman [8] and Babbs and Webber [9]. See also 
in the sections below. See [3] and [2] for approximation methods as well for 
models and processes. 

2-3 VAR 

The continuous VAR in finance was used by Lantegieg in his paper [87]. This 
process can be defined as follows. 

Definition 1.1 (Continuous- time VAR) There is a vector of state vari- 
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ables, v, that follows the process 

dv = (b + Av)dt + Gdw (1.2) 

where v is an n by 1 vector of variables, b is an n by 1 vector of parameters, 
A is n by n, G is n by k and dw is k by 1. The vector dw is a vector of 
Wiener processes, with mean 0 and variance dt. 4 

This is a slightly more general form than Langetieg. 

In Langetieg's model, the instantaneous short term rate is defined by 

r = a + p'v (1.3) 

Where (3 is an n by 1 vector of parameters and the prime indicates the 
transpose. In this notation 

n-l 

(3'v = Y,WW\ (1-4) 
Note we could also index from i=l to n instead. 

Definition 1.2 (Discrete Time VAR) The n by 1 vector y follows the 
process 

y[t] - y[t - 1] = (fio[t] + K[t}y[t - l])At + E[t] VSt (1.5) 
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The time interval can vary with time t. It is possible to allow /i,K, S to 
depend on t or be independent of t. A 



A discrete time VAR need not have an equal time step. In both a continuous 
and discrete time VAR, the parameter vectors and matrices can be indexed 
by time. Note that a continuous time VAR determines many discrete time 
VAR's, corresponding to different sets of time intervals. The matrices in the 
discrete time VAR and the probability density function, characteristic func- 
tion, Green's functions, can be obtained using the methods in the references. 

Definition 1.3 (VAR) A vector autoregression, VAR, is either a discrete 
time VAR or a continuous time VAR. 4 

2-4 Quadratic 

The quadratic VAR is the same as Langetieg's model except that r is a 
quadratic function of the state vector. 



r = a + (3'v + v'Qv 



(1.6) 
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Here Q is n by n. The expression v'Qv is a quadratic form: 

v'Qv = Y^v[^v\j\Q[i\\j} ■ (1.7) 

where the sums are from 0 to n-1 or 1 to n depending on the indexing. We 
normally use 0 to n-1. 

See Longstaff [90], Beaglehole and Tenney [19] [20], Constantinides 
Eterovic [51], Mansitre [94] Ahn, Dittmar, and Gallant [1] Leippold and Wu 
[W>], Andrew Ang and Monika Piazzesi [107] [6] (and the other papers found 
by searching on their names) and the references in the papers cited herein 
such as given in Dai and Singleton [42], [41], as well as the many found by 
searching on "quadratic term structure model" and "quadratic interest rate 
model". 

2-5 Affine and Regime Switching and Other Models 

Dai and Singleton [42], [11] review interest rate models including the affine 
models as well as regime switching models of the affine and quadratic type. 
See the text by Duffie [15]. For credit models see Duffie and Singleton [Is] 
as well as the Dai and Singleton references. See also the references cited 
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therein, and also for the affine model Merton [96], Vasicek [J. 25], Cox, Inger- 
soll and Ross [36], [37], [38], Hull and White [74], Turnbull and Milne [124], 
Jamshidian [77], [79], [80], [78] [83], [8!] Longstaff [90], Langetieg [87], Bea- 
glehole and Tenney [19], [34], Longstaff and Schwartz [92], Eterovic [51], Das 
[43] Lin Chen [32], Balduzzi, Das, Foresi, and Sundaram [11] [12], [Mi], [iu], 
Duffie, Pan, and Singleton [17], Pan [KM], and the references cited therein. 
See Shiryaev [ill] 

See Zhou [127], [J 5], Bansal and Zhou [11], Bansal, Tauchen and Zhou 
[ i 3] , and Bollerslev and Zhou [2 I] . See the work of Heath, Jarrow and Morton 
[70], [72] and models built on their approach. See Fitton [52] and Chacko 
and Das [30] for solution methods. 

2-6 Actuarial Work 

American Academy of Actuaries C3 Phase I and C3 Phase II use models for 
simulatin of assets and liabilities for determining required capital. 

See the work of Geoff Hancock of William Mercer and in Canadian In- 
stitute of Actuaries, Society of Actuaries and other meetings reporting on 
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C3 Phase II and the Canadian OSFI regulations in this area. See Geoffrey 
Hancock on bring risk into capital management [64], as well as the other 
talks in that session, all collected into a pdf file by the SOA. See [1 03], [102] 
and at SOA [I 16], [88]. 

See also the reports of C3 Phase II of the American Academy of Actuaries 
[111], [113], and [115]. See also previous reports on Universal Valuation 
System (UVS) of the American Academy of Actuaries and on Equity Indexed 
Annuities of the same organization. See the work of Mary Hardy such as 
[05] for the above organizations and in her book on guarantees in insurance 
products [()()]. A regime switching model for equity volatility and return is 
reported by Hardy and also by Hancock in the above and has been used by 
these organizations. The above reports can be found by searching on their 
titles or authors or on the web pages of those organizations. The above are 
a selection of the material of these organizations and individuals that will be 
found from these searches. 

See also Casualty Actuary Society papers and work on Dynamic Financial 
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Analysis (DFA), as well as on interest rate models and stochastic simulation. 
See also papers of other actuarial organizations including the Institute and 
Faculty of Actuaries, and the International Actuarial Association. Also the 
ETH insurance, finance and mathematics group, and the actuarial depart- 
ment at the University of Waterloo. 



2-7 Matrix methods 

See Langetieg [<S7], Beaglehole and Tenney [19], Eterovic [51], Lin Chen [32], 
[:;:;], [Hi], Duffie, Pan, and Singleton [17], Pan [10 L] Ahn, Dittmar, and Gal- 
lant [1] Leippold and Wu [«0], [fi] and the other references cited herein as 
examples of how to use matrix methods to solve for the probability density 
function over finite time intervals. See also the text by Arnold, Stochastic 
Differential Equations: Theory and Applications [7]. The book Numeri- 
cal Linear Algebra Techniques for Systems and Control by Patel, Laub and 
Doren [Ki(>] contains matrix algorithms that can be used as part of this. 
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CHAPTER 1. BACKGROUND OF THE INVENTION 
MFC VAR-ESG 



The Mathematical Finance Company VAR-Economic Scenario Generator 
(VAR-ESG) is described briefly here. This product has been used to generate 
scenarios of yield curves and equity indices. Various papers on this system 
have appeared over the years at the Society of Actuaries and Canadian In- 
stitute of Actuaries. Information is available in the SOA Library and can be t 
searched on the web for either of these. See Tenney [] I S]. Two presentations 
at conferences were [1 19] and [120]. See the important notes in these on the 
relation of this work to the work of Marjorie and Michael Hogan in the 1980's 
which was never described in any publication. The first was considered too 
difficult to print by the SOA and the second was distributed to the partici- 
pants, but not otherwise published. That book Duffie and Tavella [i')] has 
been cited in US Patent 6,173,276 by an unrelated party. Duffie and Tavella ' 
were the conference co-chairs. See the references cited therein including the 
work of Marjorie and Micheal Hogan [7^]. See Groover and Tenney [V7], 
[ r >0] Craighead and Tenney [10] The papers by Craighead and Tenney 
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and by Bobo, Tenney, and Tenney [23] are examples of application or param- 
eter estimation or development. See also papers by Tenney, Craighead in the 
SOA Library and publications, Actuarial Research Conference and ARCH, 
Canadian Institute of Actuaries, and references to them in these and other 
actuarial organizations such as the Casualty Actuarial Society. The latter 
also has extensive materials on Dynamic Financial Analysis (DFA). The SOA 
tends to call that Dynamic Financial Conditions Analysis (DFCA). 

The system in its normal mode uses the Double Mean Reverting Process™ 
(DMRP), defined as a 2 variable VAR where the interest rate is an exponen- 
tial function of the first state variable. 



with the correlation of dz x and dz 2 being p. There are different parameters 
in the risk neutral and objective measures, i.e. the Q and P measures. Here 



du — k(9 — u)dt + adz\ 



(1.8) 



dO = k 2 (02 - u)dt + a 2 dz 2 



(1.9) 
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the short rate is r = e u . Note one can introduce a time-dependent factor, so 
that r = e u a(t). 

Matrix methods referenced above are used to compute the matrices and 
vectors over discrete time intervals for calculating probability density func- 
tions and generating random variables. 

Individual yields can be modeled through solving the DMRP. One can 
also have processes on "yield residuals" such as an AR(1) process. These can 
be fit to the initial yield curve and decay according to such a process, or a 
straight line decay over a finite time interval. 

2-9 QRMC 

Quasi- Random Monte Carlo also sometimes called Low Disrepancy Sequences 
have been used with the MFC ESG based on a VAR already. This is with a 
system by Columbia University called FinDer(TM) based on the Columbia 
University patents on QRMC/LDS. See US Patent 5,940,810 Traub, Paskov, 
and Vanderhoof, and US Patent 6,058,377 Traub, Paskov, Vanderhoof, and 
Papageorgiou. 



1.3. FINANCIAL TRANSACTIONS 31 
QRMC/LDS sequences can be used to generate random numbers for both 
the regime switching and the VAR generation. These can be combined or 
separate. That is, one can form the variables to be used into one set of 
variables and apply the QRMC/LDS algorithms to that or apply it to sep- 
arate sets, e.g. the pure regime switching and the VAR. See the material 
on Columbia Universities website as well as other references such as those of 
the University of Waterloo. See the references in the paper by Tenney [121] 
on QRMC/LDS including web references. This was posted on the Canadian 
Institute of Actuaries (CIA) website sometime in 2003 and can be found by 
searching on the title or the author and several other names or within the 
organization's site. 

1.3 Financial transactions 

Financial transactions include the purchase or sale of financial contracts or 
the election of an option under such a contract or a financial option of a 
non-financial contract. They also include determining a quantity such as 
a dividend based on other financial transactions or a report on financial 
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condition. 

1.4 Basel 

The Bank for International Settlements (BIS) also often called Basel or the 
Basel Committee is an international body in the area of banking, capital 
management, risk, regulation and reporting. Its publications including its 
web site contains extensive documentation on these issues. 
Among these see 

1. Overview of the New Basel Capital Accord April 2003 

2. The Third Consultative Paper of the New Basel Capital Accord, (E), 
April 2003 

3. Basel Committee Publications in a numbered series from 1 to 106 as of 
January 2004. 

4. Basel Committee Working Papers. 

5. World Bank Basel related publications. 
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6. US Federal Reserve Basel related publications. 

1.5 OSFI 

Canada's OSFI has taken a lead in regulating financial service capital using 
simulation. 

See the memorandum headed as follows: 

Dur number: P2218-1 December 23, 1999 TO: Chief Executive 
Officers Federally Regulated Life Insurance Companies and Fra- 
ternal Benefit Societies Subject: Minimum Continuing Capital 
and Surplus Requirements (MCCSR). 

See documents of the Canadian Insitute of Actuaries (CIA) on this sub- 
ject. 

1.6 International Accounting Standards (IAS) 

See the work of the International Accounting Standards Board (IASB) espe- 
cially as it relates to insurance and fair value. 



34 CHAPTER 1. BACKGROUND OF THE INVENTION 

1.7 Financial Accounting Standards Board (FASB) 

Work on similar topics to the IAS is available from the FASB. 

1.8 Fair Value of Financial Contracts 

Fair value of financial contracts is a subject that took on a new dimension 
with the developments in modern accounting and finance. Of particular 
importance is the paper by Black Scholes [22], Garman [M], Richard [MM], 
Harrison and Kreps [G7], and Harrison and Pliska [08], [<)<)]. This also applies 
as well to papers tracing back to the portfolio methods of Hakansson's thesis - 
[!>]. These include applications to equilibrium pricing by Merton [<>>], [>'7], 
Lucas [{):}], Vasicek [125], Cox, Ingesoll and Ross [37], [38] and others. 

1.9 US House 

See hearings related information from The U.S. House Financial Services 
Committee relating to the Basel Accord on June 16, 2003, June 19, 2003 and 
other days. See also other documents on their web site found by searching 

•< 

restricted to the restriction of the .gov domain to the Senate. a 



I JO. US SENATE 

1.10 US Senate 
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See hearings related information from The U.S. Senate Committee on Bank- 
ing, Housing, and Urban Affairs "A Review of the New Basel Capital Accord" 
June, 18 2003. See also other documents on their web site found by searching 
restricted to the restriction of the .gov domain to the Senate. 

1.11 Additional Review 
11-1 Monte Carlo in Finance 

Monte Carlo was introduced into finance for the valuation of securities by 
Boyle in 1977 [27]. Since then it has been widely used for both valuation and 
risk analysis. 

11-2 Option Pricing 

Option pricing formulas were developed by Bachelier [10], Kruizenga [S<>], 
Sprenkle [112] and Boness [2. r .] J . Black-Scholes [22] showed how to derive 

a formula of the same form as Boness [2- r >] using an equilibrium approach. 

1 Some notes in the next few subsections are based on those in Davlin et al. Tenney 
[ i] An earlier version of this paper was posted on the SOA website. 
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Boness and the others had been unable to come up with the right discount 
rate and expected rate of return, the risk free rate in both cases, these were 
discovered and proven by Black-Scholes . Black-Scholes report a no-arbitrage 
derivation based on a suggestion by Merton who showed the Boness type 
formula with the Black-Scholes parameter restrictions for equity options also 
obtains if interest rates are normally distributed. Black-Scholes showed how 
to derive a partial differential equation of the McKean-Samuelson type [H-j] 
[109]. Merton [98] extended that equation to include a second random source 
from interest rates for pricing equity options only. Cox-Ross [39] showed how 
to interpret the Black-Scholes solution in terms of risk neutral probability 
and to price options with other stock price processes. This was an economic 
interpretation of the mathematics of Black-Scholes and Merton which was 
already risk neutral probability. 

The general approach to arbitrage in a single currency was done by Gar- 
man [r>;j] at the same time as Richard developed it for just random interest 

rates [108]. 2 The Garman approach was for any type of security contingent 
2 This was in a 1976 Carnegie Mellon working paper. 
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on any type of random variable. 

11-3 Interest Rate Models 

Monte Carlo is the leading method used to price mortgages and equilibrium 
mortgage interest rates are determined by the use of Monte Carlo simula- 
tion models. The mortgage applications use interest rate models as the main 
random factor. One factor interest rate models under equilibrium were de- 
veloped by Vasicek [125] and Cox, Ingersoll and Ross [38] 3 . Richard [108] 
and Cox, Ingersoll and Ross [38] extended the one factor CIR model to two 
factors independently. 

The key base of modern multi-factor interest rate models with closed form 
solutions are the multivariate normal models of Langetieg [b7] and the two 
factor square root model of Richard [108] and Cox, Ingersoll and Ross (CIR) 

M ■ 

Special cases of Langetieg's model are Ho-Lee, , [72], Jamshidian, [77] 
[7!)] [SO] , and Hull and White and [71] . 
3 Both were working papers in 1976 
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Later models built on Langetieg include the quadratic model Longstaff 
[HO], Beaglehole et al, [19] [20] , Jamshidian [81] 4 Eterovic [51], Ahn, Dittmar 
, and Gallant [1], Leippold and Wu [89] and others, Lin Chen's [:V2] 3 factor 
model and the affine model of Duffie and Kan [16], and the Heath, Jarrow, 
and Morton (HJM) methodology [70]. 

Langetieg derived a Boness-like formula for options on stocks with his 
interest rate model applying Merton's argument on Black-Scholes. Hull and 
White showed that Langetieg's approach could be slightly modified to ap- 
ply to their version of the Langetieg model. State prices or option prices 
for these models were developed by Cox, Ingersoll and Ross [3N] Jamshidian - 
[79] , Longstaff [9.1], Hull and White [74], Beaglehole et al. [19], Milne and 
Turnbull [TJ-1], Chen and Scott [M l], Longstaff and Schwartz [0?], Constan- 
tinides [:W], Lin Chen [;i2], Manistre [91], Duffie, Pan and Singleton [i7] and 
others. Earlier work on state prices or Green's functions in finance traces 
back to McKean [05], Garman [!">;>], Cox, Ingersoll and Ross [■>**], Breeden 

and Litzenberger [29], Banz and Miller [in], Ingersoll [75], and Merton [<»■)]. 
4 See Manistre pM] for a general discussion of multifactor models. 
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In addition, Lucas [93] derived state prices building on the foundation of 
Hakansson's [G^], [oiJ] approach to optimal savings, consumption and portfo- 
lio choice. Lucas developed a state price using marginal rates of substitution 
that arise in the Hakansson type optimizations of consumption, savings and 
portfolio decisions in a multiple-period context. Cox, Ingersoll and Ross and 
Vasicek's equilibriums are also built on the Hakansson foundation. 

Because of the prepayment models and the frequent use of interest rate 
models without closed form solution, Monte Carlo and Low Discrepancy 
Sequences are often used for mortgage pricing and risk analysis. An example 
of a model without a closed form is the DMRPTM 5 model, see for example 
Groover [ r n], Groover et al. [ r ;7] and Craighead et al. [10]. 

11-4 Insurance Applications 

Monte Carlo has been applied to insurance products or asset liability man- 
agement by Boyle, Brender, Brown, Craighead, Embrechts, Fitton, Hancock, 

, Hardy, Manistre, and Panjer in numerous papers and actuarial studies to 
5 Tradernark Mark Termey and Mathematical Finance Company 
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name but a few. It has been used to set Segregated Fund Guarantee reg- 
ulations in Canada. It has been used for interpreting CARVM for Equity 
Indexed Annuities, for C-3 Phase One risk based capital for interest rates, 
C-3 Phase Two for variable annuity guarantees, and other studied by the 
American Academy of Actuaries and Canadian Institute of Actuaries . 

Monte Carlo has been applied to guarantees on segregated funds by Mary 
Hardy using a regime switching approach [65]. Error bounds for this model 
or related ones can be calculated with discrepancy. Low Discrepancy Meth- 
ods or Quasi-Random Monte Carlo can be combined with that model for 
analyzing insurance problems such as variable annuity guarantees including 
reserve and capital requirements. 

11-5 Monte Carlo 

Monte Carlo was developed at Los Alamos in the 1940's for calculations re- 
lated to nuclear physics. It has since been applied to a variety of problems 
from analyzing scattering experiments in accelerators to financial applica- 
tions. 
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11-6 Quasi- Random Monte Carlo 

Quasi-Random Monte Carlo means deterministic sequences or points. They 
are usually chosen because they have a lower discrepancy than "randomly" 
chosen points. Paskov and Papageorgiou [105] review the history of apply- 
ing Quasi-Random Monte Carlo (QRMC) to high-dimensional problems in 
finance. The text by Niederreiter [101] is one of the classic texts on QRMC 
and is the basis of most of the theorems presented later. 

Some recent applications of QRMC have been made by Boyle, Broadie 
and Glasserman [28] and Boyle and Tan [2(>]. 

1.12 Copulas 

See in particular Chapter 2 of Nelsen [100] and Chapter 2 of Embrechts, 
Lindskog and McNeil [f>0]. The errata for Nelsen's book is at his university's 
webpage. (search on the title or author). However, the entire documents are 
important material on copulas. A good introduction to applying copulas to 
reinsurance is by Gary Venter [ I *JI>], This has many good pictures of copulas. 
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This is available at the casualty actuary website. 

1.13 FX modelling 

FX stands for foreign exchange. See Davlin and Tenney [14]. An earlier 
version was posted on the SOA library website. The following references are 
taken from that paper. 

Krugman The Garman-Kohlhagen model without stochas- 
tic interest rates appeared in 1983 [51] and also derived a Boness 
type formula for options on exchange rates. The arbitrage free 
framework for pricing pricing with stochastic interest rates and 
exchange rates appears in Ingersoll [7G], Amin and Jarrow [I], 
Jamshidian [8:2]. The material of section two was based on some 
fragmentary notes [117] that appeared in a fragmentary draft pa- 
per by Beaglehole et al [17]. Some related work on stochastic 
exchange rates and stochastic volatility are the notes of a paper 
on stochastic volatility by Beaglehole et al. [J..s] 5 Heston [71] and 
Knoch [SI], Melino and Turnbull [122] and [122,]. 
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13-1 MFC-FX 

A system was used by MFC where a variable v was constructed as follows. 

v = 7 + a £(Ar«)At + /?w (1.10) 
t 

where Ar t is the difference in interest rates in two currencies at an observi- 
ation point and w is some element from the VAR. These were obtained by 
exponentiating elements of the VAR in each currency. The v variable is con- 
structed from some reference point in time, or a fixed lag. Let Af be the 
change in the log of the exchange rate. It is modelled that 

Af = a + bv (1.11) 

A continuous time version of this is easily constructed. 

dv = <pdt + a(Ar t )dt + pdw ' (1.1-2) 

(in the case that a, (5 don't vary with time). Depending on the specification 
of the exchange rate and interest rate differential, an Ito type term may have 
to be added or subtracted or such a term as occurs in the derivations in the 
Davlin Tenney paper. 
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1.14 Financial Advisory Systems 

See U.S. Patents No. 5,918,217, 6,012,044, 6,021,397, 6,125,355 and 6,292,787; 

1.15 Other financial systems 

See U.S. Patent No. 5,193,056 on a "Data Processing System for Hub and 
Spoke Financial Services Configuration." 
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37 CFR 1.77 b "(7) Brief summary of the invention." Starts here. 
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Chapter 2 

Summary of the Invention 



A regime switching vector autoregression (RS-VAR) is defined as a vector 
autoregression in which the parameters of the vector autoregression are func- 
tions of a set of discrete indices, which consitute the regimes. This process 
can be applied to interest rate models, default models, and other financial 
models. This can be done in the "objective" or P-measure or the risk-neutral 
or Q-measure of finance or other measures. One set of applications include 
calculation of prices, cashflows, capital, reserves, defaults, and other vari- 
ables. Another set includes transactions using these including purchases and 
sales, producing and/or sending reports, advisory services, portfolio strategy, 
etc. 

47 
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Typically, these applications involve using technical means such as com- 
puters or the internet. An additional class of applications are portfolios and 
financial products made with these methods. 



Chapter 3 

Brief Description of the 
Drawings 

3.1 RS-VAR generation 

Figure 1 is a flow diagram of the generation of the Regime Switching Vector 
Auto- Regression (VAR). 

3.2 Discretized RS-VAR generation 

Fig 2 is a flow diagram of the Discretized RS-VAR generator Time Loop. 

3.3 Discrete Time Regime Probability Prop- 
agation 

Fig 3 is a flow diagram of Discrete Time Regime Probability Propagation. 

49 
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3.4 Matrix times random vector 

Fig 4 is a flow diagram of a matrix times random vector calculation in the 
discrete time portion of the VAR vector calculation. 

3.5 State Variable Fitting 

Fig 5 is a flow diagram of the State Variable Fitting. 

3.6 Discretization 

Fig 6 is a flow diagram of calculation of the mean, variance-covariance matrix, 
the standard deviations and the correlations of the discretized process. 

3.7 Yield Grid Generator 

Fig 7 is a flow diagram of the Yield/Price Grid Generation. 

3.8 Yield Scenario Generator 

Fig 8 is a flow diagram of the Yield Scenarios Generation. 
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3.9 Basic ESG 

Fig 9 is a flow diagram of the ESG. 

3.10 Auxiliary Scenario Processor 

Fig 10 is a flow diagram of an Auxiliary Scenario Processor. Examples are 
combining yield grid data and the state variables and also auxiliary data of 
some sort being tracked. 

3.11 Statistical Analyzer Module 

Fig 11 is a flow diagram of the statistical analysis of the scenario data. 

3.12 Price 

Fig 12 is a flow diagram of the price calculation using scenarios. 

3.13 Sensitivities 

Fig 13 is a flow diagram of the sensitivies calculation. This includes price 
derivatives with respect to state variables or other paramters such as a spread 
or option adjusted spread. It includes durations with respect to the state 



52 CHAPTER 3. BRIEF DESCRIPTION OF THE DRAWINGS 

variables and other parameters and a convexity matrix with respect to them 
and other parameters. 



3.14 Reserves 

Fig 14 is a flow diagram of reserves calculation. 

3.15 Capital 

Fig 15 is a flow diagram of capital calculation. 

3.16 Portfolio Simulator 

Fig 16 is a flow diagram of a portfolio simulator. The portfolio can include 
assets and liabilities. 

3.17 portfolio simulator scenario loop man- 
ager 

Fig. 17 is a flow diagram of a portfolio simulator scenario loop manager. 
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3.18 portfolio simulator scenario time loop 
manager 

Fig 18 is a flow diagram of a portfolio simulator scenario time loop manager. 

3.19 Preparation for Transaction End Use 

Fig 19 is a flow diagram of preparation for a Transaction End Use. 

3.20 Transaction End Use 

Fig 20 is a flow diagram of execution of a Transaction End Use. 
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Part I 
Detailed Description 
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Chapter 4 
Simulation 

4.1 Introduction 

In simulation, we start out by thinking we have the true model of nature or 
the economy or the economic process and we know everything. We don't have 
any observation problems, we know the parameters, we know the models, we 
know how to simulate. 

We have as our theoretical model the idea that there are hidden variables 
which we can model as coming from such joint distribution like the multi- 
variate normal. We use the latter because it has nice properties especially 
for modeling many random variables at once. 

We have a model though that individual data series that we see in the 
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economy, like the S & P 500 or the price of Ford stock are not normally 
distributed nor are simple transforms of them like the difference in natural 
logarithms of them at adjacent time points. Nor is the arithmetic return 
from one month end to the next month end normally distributed. 

So we have a theoretical model that we can simulate with which starts 
out with the hidden multivariate normal distribution of many variables and 
simulates those first. Then for each such variable we transform it so as to 
produce the observable variables such as Ford stock price or the return of Ford 
stock price which we also classify as an observable because the transformation 
involves no estimated parameters. 

So step one is generate all the normal variables at once for a given time 
period using a multivariate normal distribution. This models their correla- 
tion. Now we take each individual normal and turn it into an observable. 
The way we do this is by matching the univariate normal distribution's cu- 
mulative distribution to the univariate cumulative distribution of some other 
distribution like logistic. At this point, we know the parameters of the logistic 
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or final distribution as we might call it or final output distribution. 

We also know the mean and variance of the inner core normal. In this 
model, it doesn't actually matter what those are as long as the variance 
isn't zero. The reason is that the parameters of the "output" univariate 
distribution will undo any effect of the mean or non-zero variance of the 
normal. 

So we say that the cumulative normal (of the univariate normal distribu- 
tion) up to the simulated value x equals the cumulative distribution of the 
"output" distribution, say the logistic, whose value we call y. The distribu- 
tion of y values doesn't depend on the mean and variance of the normal as 
long as we use. the same mean and variance for x's cumulative distribution 
to transform x to y as we used to simulate x. 

We thus solve for y by equating its logistic or whatever cumulative density 
to the cumulative density of x, with the same mean and variance for the 
cumulative of x as are used to generate x. And as long as that is done and 
x is univariate normal, it doesn't matter what mean and variance were used 
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to generate x as long as the variance wasn't zero. 

So we solve for y and that is now our output variable. The variable y 
might however be something now like the change in the logarithm of the 
stock price or index over some period of time, like a month. 

Note that the transformation has to pick a specific time interval for it 
to be based on which has some specified logistic. It is possible to define a 
continuous in time output variable but it will only be logistic with specified 
parameters over some specific interval. 

4.2 MVN 

Let x be an n by 1 column vector. The 1 by n row vector x' is given by 

x' = (x u ...,x n ) (4.i) 

We assume that x is multivariate normal, MVN, with probability density 
function 
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4.3 Univariate Distributions 



For each x { we transform to a variable y u defined by the following procedure. 
Let jii and a* be the parameters of the univariate distribution for X{. Let 
Fi(xi] fi, a 2 ) be the normal cumulative distribution with mean fi and variance 
a 2 for variable X{. We have 



Let Gi(yi]rj) be the cumulative distribution function of a single variable 
with parameter column vector rj of dimension ki by 1. Let the corresponding 
probability density function be g{yi\ 77) , with appropriate modifications made 
for point masses, such as including a 5 function in g. Let yi be determined 
by the equation 



We solve the above equation for conditional on the parameter vector rj 
after simulating the vector x including the value Xi. In this approach, d can 
be any univariate distribution 




(4.3) 



Fi(xi;n, v 2 



') = Gi( yi ;v) 



(4.4) 
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4.4 Discrete time simulation 

We can simulate in discrete time, most simply when the time intervals are 
equal. 

4.5 Combining with other variables 

We can combine the variables generated in this manner with other variables, 
which themselves may follow a continuous time process. The vector x of 
MVN for tranforming to logistics can be generated as part of that process, 
and we can for example choose to make that subvector have a mean of zero 
and standard deviations of one over the time intervals at which we wish 
results. 1 This sub- vector then furnishes a vector of MVN for the above 
transformation. For the case of fixed time intervals for the desired logistics, 
we simply specify the corresponding logistic distribution parameters for an 
interval of that size and use the generated MVN over that period. 



This is not really necessary. 



Chapter 5 

Esimation of Parameters 



For estimating parameters we have an inverse direction from simulation. In 
the case of simulation we start first with the simulation of the multivariate 
normal vector of variables first. Then we transform those to the individual 
output variables. This can be done by a transformation using cumulative 
density functions of the normal and any univariate distribution, such as a 
logistic. 

In the case of estimation, we start with data that is already transformed 

so to speak by "nature" or the economy. Our job as experimentalists or as 

econometricians is then to reverse engineer nature. So we get from nature or 

the economy, variables that individually we identify as logistic or whatever. 

65 
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We believe there is an inner core of correlation that is hidden by the 
logistic or univariate distributions of individual data series, such as the S & 
P 500, Ford stock price, etc. 

We can't observe the correlation matrix directly nor estimate it directly 
from the observed univariate series, e.g. Ford stock price at month end, S & 
P 500 at month end from January 1950 to January 2003, or whatever period 
we use or data series we observe. 

Because these invididual data series are logistic or whatever, we can't 
directly estimate the correlation matrix of them. The multivariate logistic is 
a very constrained distribution that we don't want to be restricted by. We 
might also want to use other univariate distributions than logistic. 

So to uncover the hidden correlation of this collection of individual series, 
modeled by a collection of independent or seemingly independent univariate 
distributions, we want to transform each individual data series so that the 
transformed series is univariate normal. We then make the leap of faith, 
or assumption that these transformed series are in fact multivariate normal, 
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MVN, and that we can estimate their correlation matrix by standard means 
for a MVN given a collection of individual normals that are further assumed 
to be multivariate normal. As is well known, that is not guaranteed, but as 
modelers of nature or the economy we make that our assumption because its 
easier to build a model. If that doesn't work we can try something else. 

So as econometricians our first problem is to estimate univariate distri- 
butions for each series that comes to us from nature or the economy, the 
SPX, Ford, whatever. We in fact first transform this by using differences in 
the natural logarithms of the series or constructing the arithmetic returns. 
Because these transforms involve no parameters to be estimated we think of 
this as still being the data that comes to us from nature or the economy. 

5.1 Estimation of Parameters 

Suppose we are given a vector of data y t , t=l,...,T, where y t is a vector of 
dimension n by 1. For each i, we take the time series y ti , t=l,...,T, where 
y ti is the i-th component of y t . Given this data we can estimate a univariate 
distribution separately for each i. This can be done using standard univariate 
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methods. Let the estimated parameters for each distribution be given by the 
vector rj h where rji is a k { by 1 column vector. 

Given these parameters, we can now transform to a distribution that is ' 
univariate normal. We discuss two ways to do this. One is to specify that 
the univariate normal has mean zero and variance 1. We then solve for x it 
for t=l,...,T such that 

/ X it 1 
j z 7wr {,) =G,(B "^ <51) 

This gives us a set of vectors x u t=l,. ..,T. We now estimate the correlation 
matrix among these using standard methods. Let this correlation matrix of - 
the vector x t at any time point be denoted W \ 

We simulate with this approach a MVN with mean zero, and unit stan- 
dard deviations with correlation matrix W. We then transform each y± indi- 
vidually using Gi(.,bi). 

As an alternative, we could specify a vector fx for the MVN and a vector 
of standard deviations a. These could be the values estimated in the sample. 

As a short cut to estimate the correlation matrix, we can use the correla- 
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tion matrix of the raw data. When dealing with some data, we may wish to 

take a transformation like the natural logarithm first. If we deal with stock 
data, for example we can first take the change in the natural logs of the stock 
prices and treat those as the y's. In this form, we may use the correlation 
matrix of the y's as a quick estimate of the correlation matrix of the normal 
distribution. 
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Chapter 6 
Statistics 



6.1 Univariate Logistic Distribution 

For any y between -co and oo, let the cumulative distribution of y given 
parameters a and (5 be 

G(y,*,P)= 1 + e _\ x _ aW (6.i) 

This can also be written 

G{y\ a, 0) = |(1 + tanh(~(x - a)//?)) (6.2) 

G is the cumulative logistic distribution and y is said to be logistically dis- 
tributed or to have the logistic distribution. 
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6.2 Logistic Normal Conversion 

To do conversion we equate the cumulative distribution function of the Nor- 
mal to that of the Logistic. 



If we are simulating we simulate x with mean 0 and variance 1 and then 
use our model assumptions for a< and ft for series i, e.g. Ford or S & P 500 
or whatever. 

If we are estimating, we first estimate a, and ft from a time series of y it 
and then we convert to x using the above formula given the observed value 
of Va- 
in both cases y it is the return, log or arithmetic return, not the stock 
price or index value itself. For example x and y are both varying over -oo 



The above normal has zero mean and variance 1 which may seem counter- 
intuitive. Why doesn't it have a mean and variance corresponding to that of 
the logistic? 




(6.3) 



to oo. 
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The reason is it wouldn't matter. Suppose we use some n and a in the 



ulation as well as conversion. As long as we simulate with the same fi and 
a that we convert with, it doesn't matter what n and a are as long as a is 
non-zero. We could use the values estimated from data, but we can also just 
use fi = 0 and a = 1. 
2-1 Proof 

Suppose that we started with price data, then converted to differences of 
logs and estimated a logistic on that data. Suppose we also estimated the 
mean and standard deviation, and those had values /x and a. We have the 
2/it,t=2,..,T, and we solve for x it , t=2,..,T. 



transformation 




(6.4) 



won't make any difference once we look at the combined effect of si 




(6.5) 
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We then estimate our varianee-covariance matrix W. However, for now 

lets pretend we have just one series i. Now we come to the simulation point. 

And imagine that we simulate u it using a normal distribution with // = 0 

and (7 = 1. 

Now we could calculate 

x it = n + auit (6.6) 
We now calculate the cumulative distibution function, F using 

dz . e-fo-dV" 3 ( 6 .7) 

We can substitute for x it as 



. iz jm?*-^ (68) 



Now to calculate this integral, we tranform z to a standard normal. We do 
this by setting 

v = (z~ h)/<j (6.9) 
Making this transformation, we find that 

F °fy^r" (6io » 
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Now this is the same thing as if we had just simulated u it and then used this 
cumulative distribution to calculate F. 

The value of yu depends on F, a* and $ not on how F was calculated, 
as long as it produces the same F. So if we simulate u as mean 0 variance 
1, and then transform u to some x which has mean /i and a and calculate F 
using the same /j and a its the same as if we calculate F using u and using 
// = 0 and cr — 1. 

2-2 Extension to other distributions 

The proof just given that the conversion from the normal to the non-normal 
distribution was independent of the normal's mean and standard deviation 
as long as the same parameters of mean and standard deviation are used to 
generate the normal deviates as to do the transformation did not rely on the 
characteristics of the logistic, and the same proof applies to any non-normal 
distribution. 
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6.3 Gamma Distribution 

Reference Morris DeGroot Probability and Statistics p236. The probability 
density for the Gamma Distribution is 



g(x\aj) = -^»-V* (6.11) 



for x > 0 and 0 otherwise. Here T(a) is the Gamma Function, which gener- 
alizes the factorial, T(n) = (n — 1)!. 

This gives us a 2 parameter distribution. We can compute the cumulative 
distribution function as 

G(y\a,0)= [ V dxf(x\a.J) (6.12) 

JQ 

So that 

G(y\a,0)= I" dx^-x^e-^ (6.13) 
Jo r(a) 

We then convert from normals to this Gamma by first simlating a normal 
deviate u with cumlative density 

F{u-^) = f U dz-=±=e-^-^ 2 l^ (6.14) 
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and then converting to y such that 

G(y\a 1 0) = F{u ]f j, ) a 2 ) (6.15) 

We have to interpret the results somewhat differently though. Since y is 
between 0 and oo it can be interpreted as the new price divided by the old 
price. In the case of the Logistic Distribution, we interpreted y as the change 
in the log of the stock price. Now we have to interpret it as the new stock 
price divided by the old stock price. 

6.4 Beta Distribution 

The beta distribution has probability density function 

(616) 

for 1 > x > 0 and 0 otherwise. 

To represent stock returns with this distribution it is necessary to convert 
the stock return or stock price to a variable between 0 and 1 and vice versa. 
We can do this by taking a lognormal variable and converting it into a unit 
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deviate by calculating u from x by 

u = F(x;fi T ,a T ) (6.17) 

or given x, solving for u. The superscript T reminds us these are the trans- 
formation parameters. One choice is fi T = 0 and a T = 1. 

Here F is the cumulative distribution of some lognormal variable. This 
can differ from the normal used for the generation of the random normals. 

If we generated a normal x, converted that to u and then from that u to 
y for the Beta and then from that y back to x' as a normal, then we would 
have that x'wasa normal determined essentially by the initial x directly. So 
we need a different method to generate a distribution of data. 

4-1 Alternative generation of multivariate vector 

We can instead use the univariate distributions of logistic, beta, gamma, etc 
and generate uniform deviates u { . We can then transform these to normals 
x { . We can then take a linear combination of these x u yj = £ R jiXi , where 
R is some matrix. Note the dimension of the x vector can be larger than the 
dimension of the sub- vector of y we wish to use. We can now treat the y s as 
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changes in logs of stock returns. Or we could transform them to be between 
0 and oo and interpret those as ratios of prices. 

In this approach we have a fundamentally discrete time generation pro- 
cess. 
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Chapter 7 
Regime Switching 

7.1 2 Regime Switching Model Lognornmal 



For this model, let the parameters be fa and Oi in state i, i=l,2. The param- 
eters of the model shift at random times. Let p iS be the transition rate per 
unit of time of going from state j to state i. 

The lognormal model is 

dx = Hidt + Gidz (7.1) 

when in state i. 
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7.2 Estimating the parameters 

The parameters are estimated using maximum likelihood as in Hamilton, 
Time Series Analysis, Chapter 22. See for example Mary Hardy's North 
Ameican Actuarial Journal paper for an application to the two state regime 
switching lognormal. 

7.3 Estimating initial state 

The initial state can be represented as known in one of the regimes, or as a 
probability density over them. 

3-1 maximum likelihood 

Can use the maximum likelihood method in Hamilton. This produces the 
probability density. 

3-2 Trailing Data 

We can also use trailing data, especially daily data if available. We estimate 
the parameters in the previous N days, and then use a simple comparison to 
the two states of the model and use the one that is closest. Alternatively, a 
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weight can be given to each of the possible states, e.g. if the trailing vol is 
20 percent and there are two vol states of 15 and 25 percent, one could use 
a weighted average, computed in standard deviation or variance for example 
as the measure. 



7.4 Regime Switching Non-normal Models 

We can generalize the discussion of the last section to cover the other uni- 
variate models besides the lognormal. Let the model have a parameter vector 
77, with elements ^, i=l,--,n ??J where are the number of elements in the 
vector 77. Let there be n r regimes of the parameters or discrete parameter 
sets. 

Using the same methods as for the lognormal model other models can be 
estimated. This is done using maximum likelihood as in Hamilton. 
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7.5 Advantages of Regime Switching 



5-1 Data description 

There Regime Switching Model (RSM) describes the effect in data of periods 
of different values of parameters. Periods of high volatility in particular are 
important in modeling indices, especial equity indices. The Regime Switch- 
ing Model does more than model the -fat tails of the data, it also models the 
tendency for one fat tail event to be followed by another, i.e. a switch in 
regime in which large losses or gains are more likely. This is important to 
option writers because such gains or losses have a non-linear effect on the op- 
tion price. This non-linear effect and the probability of greater price changes 
then results in diffrent prices for the options, so that their expected return 
in the risk neutral probability is the same. In particular, the optionality of 
the option, its convexity is worth more when volatility is higher. This leads 

* 

to a higher price of the option, put or call. 
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5-2 Risks to Option Writers Modeled by Regime Switch- 
ing 

Option writers, such as by variable annuity sellers, have the risk that volatil- 
ity will change from what it is or they expect. When this happens, the Greek's 
of the options written change. This means that after a change in parameter 
values, the change in option prices for a given change in the stock price will 
be different. As a consequence, hedges set up to match Greeks, i.e. to match 
risk of the asset portfolio to the liability portfolio will no longer match. This 
is because the. asset portfolio was chosen for a given set of Greeks. When the 
parameters change the asset and liability portfolio options all change in their 
Greeks, but because the Greeks are non-linear functions of the parameters 
and the mix of options in the asset and liability portfolio are different the 
options will change their prices from the parameter change and their Greeks, 
so that there is a discrete change in price from the parameter change and a 
discrete change in Greeks. This can cause a loss itself and also means the 
asset portfolio no longer hedges the liability portfolio. 
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The changes in prices and Greeks from the parameter change is a risk that 
can be analyzed using theregime switching model. Moreover, the change in 
hedge portfolio after the jump to one more appropriate to the new regime of 
parameters can be analyzed using simulations based on the regime switching ' 
model. One can estimate the cost for example of changing the hedges after 
a parameter change. Moreover, one can estimate the cost of the parameter 
changes and develop a portfolio to partially hedge this. 

7.6 Canadian and US Insurance Regulation 

Canadian and draft US C3 Phase II regulations for capital have a simulation 
option or requirement (use of tables is allowed in Canada). This does not - 
require use of regime switching, although the calibration table in Canada 
was developed using regime switching analysis. Distributions with fat tails 
meet the calibration table with parameters that don't have super fat tails. 
By using a high enough volatility the lognormal will meet the calibration 
table. The calibration table requires the probability for large returns or large 
losses of specified size to exceed certain minimum probabilities. If one uses a 
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lognormal and increases its volatility to fit the table, to meet all the points 
one needs a very large volatility, so that one exceeds the required probability 
for several of the points in the table. With the regime switching or other fat 
tailed distributions one can come closer to just matching each point in the 
calibration table. 
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Chapter 8 
ESG 

8.1 DMRP 

Let u and 9 follow the process 



du = ki(9 - u)dt + aidzidO = k 2 (9 2 - 9)dt + cr 2 dz 2 (8.1) 

under realistic probability. The correlation between dz x and dz 2 is p. Let the 
risk neutral version be 

du = 4(0 - + <j x dzid(j) = 4(02 - + a 2 dz 2 (8.2) 

Where the variable 0 is related to 0 by a linear transformation. The correla- 
tion between dz x and dz 2 is still p and the values of o\ and a 2 are the same 

in real as risk neutral as is required from no-arbitrage. 
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The instantaneous short term interest rate is 

r = e u (8.3) 
The prices of bonds are solved from 
)f\B uu +B u K\{<j>-u)+^^ = 0 (8.4) 

The parameters of the real and risk neutral processes are estimated from 
historical data over different time periods starting in the 1950's to date for 
U.S. Parameters for other countries and for the Euro are also available as well 
as dual linked currencies such as US Yen or Euro Yen in which the exchange 
rate is part of the state vector v and multiple sets of u, 6 or u } <f> are used for 
each economy. 

In simulating real probability scenarios we can simulate u and 0 but with 
realistic probability parameters for this process, yet a 3rd process from the 
above. Given u and 0 we then calculate bond prices and yields. Alternatively, 
we can simulate u and 6 and then calculate 0 and then calculate bond prices. 
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8.2 ESG™ Simulation 

We first consider the case without regime switching. Let the vector v be m 
by 1, and let its elements be v 0 == u, v 1 = 0, let v i+1 = x u for i = l,..,n where ■ 
n is the dimension of the x vector, and we index x starting from 1. We now 
form the process 

dv = (b + Av)dt + Gdw (8.5) 

where b is m by 1, A is m by m and G is m by m. The vector dw is m by 
1 and is a vector of mean 0, variance dt independent Wiener processes. We 
simulate v over some set of time intervals, which we suppose for simplicity 
are of equal size, At. We assume we start at t = 0, and generate a time series 
v u t=0,..,T. So we simulate v Xj ...,v T and we start with v 0 . For stock return 
elements, i > 1 % = 0 and for the interets rate model we use v 0Q = u and 
Vqi = 9. Here the first index of v is the time index t, and the second index 
the index j for the vector component. 

The first two elements of b and the first 2 by 2 sub-matrix of A are formed 
so as to replicate the DMRP process above. All other elements of A and b 



92 CHAPTER 8. ESG 

are zero. 

The correlation matrix GG is estimated for the joint correlation matrix 
of u, 6 and x, i.e. v. 

We now simulate v and for each element of x, we have x { = v^ x suppress- 
ing the t index, and remembering we are indexing x starting with 1 and v 
starting with 0. 

We now convert each x; to y 4 using the univariate distribution for This 
can be normal or logistic or some other univariate distribution. We use the 
appropriate parameter vector b { to do this, taking care our parameters are 
approriate for the time interval A. As proven earlier it doesn't matter what 
the mean and variance of the x { are at each t, as long as we use the same 
ones for the conversion as for the simulation. We can in fact just use unit 
normals used in generating the w[t] vector. 

8.3 ESG with Regime Switching 

We have the vector of continuous variables 

dv = (b + Av)dt + Gdw (8.6) 
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In addition to this is the state index s, which we interpret as a set of a finite 
number of states. We have a mapping from s to b,A and G, assigning for 
each s, b(s), A(s) and G(s). Different values of s can then correspond to 
the different paramters of the univariate models for each index. We have a 
process on s, of the form p(s\s)dt where p{s\s) is 4 the transition rate per 
unit time to go from state s to state s\ In addition the transformation 
functions from v to y are indexed by s in addition to their other parameters 
so y = h(v,s,t) is the transformation function from v to y contingent on 
state s. 

8.4 Indices and Interest Rates 
4-1 Risk Neutral 

Suppose we have a stock index with cumulative total return value V. So at 
time t=0, we have V=0. 

Consider the case that the index's returns are lognormally distributed. 
The risk neutral process on V is 

dV = rVdt + crVdz (8.7) 
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where dz is correlated to the other elements of the vector y. 

Over a finite interval from t to t+h, we have 

l nV [t + h}- lnV[t] = J[r- l -a 2 \di! + J adz (8.8) 

lnV[t + h)- lnV[t] = J r{t')dt' - ±a 2 h + J adz (8.9) 
where we assume that a is constant over the time interval h. 
If r(t') was constant, then we could do the integral on r and get 

lnV[t + h]- lnV[t] = r(t)h - ^a 2 h + J adz (8.10) 
However, for the DMRP model, r(t) is not constant. Instead it is a 
random process over the time interval. At time t, we want the prospective 
distribution of V[t+h] at the end of the time interval [t,t+h] conditional on 
information at time t. We can approximate the above process in risk neutral 
by 



lnV[t + h}- lnV[t] = C - \a 2 h + ! adz 



(8.11) 
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where C is the yield on a bond from t to t+h, conditional on the state 
vector y at t. So 



C -C fc = B{t,t + h,y) (8.12) 

where B(t,t+h,y) is the zero coupon bond price at time t that matures at 
time t+h conditional on the state vector y at time t. 

4-2 Realistic Probability 

We can write the return on the index in realistic probability as 

. lnV[t + h}- lnV[t] = a + bQ - ^a 2 h + J adz (8.13) 

Where a and b are constants estimated from a regression or by theory. 
For example, one could use the CAPM values. This could be represented as 

^ = ( + ^1-0 (8.14) 

where (3 is the CAPM beta of the index relative to some market index and 
( is the yield over the interval h on a zero coupon bond, and E[R m ) is the 
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expected return on the market index. 



So for the CAPM one chooses 



(8.15) 



6=1-/3 



(8.16) 



4- 3 Bond Index 

The relation of bond indices to interest rates is discussed in a separate section. 

8.5 ESG Estimation 

5- 1 Correlation Matrix 

To estimate the correlation matrix of the ESG, we proceed along the lines of 
Chapter 2, but now include past values of u and 6 estimated from the yield 
curve on past dates to get these values and then estimate the correlation 
matrix of u, 9 and x where x is the state vector y not including u and 9. 

As before, we estimate each time-series in x, e.g. S & P 500 or Ford 
stock price first as a univariate distribution, e.g. Logistic and then use the 
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estimated parameters to transform to the normals using the cumulative dis- 
tribution method. 

Suppose we have a set of indices, i=l,...,JV>. Let the total return on each 
index over a set of time intervals be v(i,t k ) for index i at time t k . We use 
the total return from tk-i to t k . 

J') = - KOIbO'.*) - W))/D (8.17) 

t 

D is an appropriate divisor, such as N k or N k - p where p is the number of 
parameters estimated and N k is the number of time points for which there 
are total returns. Note if we start with index values this is one less than 
the number of time points for which we have index values. Here v(jj is the 
average of v(i, t k ), k = 1, N k . 

8.6 Index Simulation 

A set of Nj indices are selected. The state vector y is now N y = 2 + Nj in 
dimension, unless there are additional elements in it. 
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6-1 Input 

The current value of the index is usually input as 1 instead of its actual ' 
numerical value. 

6 — 2 Process 

The change in the index is lognormal in the simplest case. We find the 
change in the log of the index as the sum of its expected change and its 
random change. These come out of the Vector Autoregression (VAR) model 

* 

on the vector y. 

dy = (b + Ay)dt + Gdw (8.18) 

Here b,y and dw are vectors of dimension N y by 1. A and G are matrices of 
dimension N y by N y . 

Each of the index returns is lognormal conditional on the starting value 
of y. We can however, use the regime switching model to augment this. 
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6- 3 Output 

The output is the log of the index value or its transformed value. 

8.7 Estimation of Index Returns 

Index returns are estimated using past data. A set of index returns over a 
period of time are converted into log returns by transformation. These are 
then correlated to each other and to the u and 8 variables of the model. The 
standard* assumption is to assume that there is no dependence on y in this 
return, i.e. the A matrix has zeros in the row corresponding to the index. 

SPX data at monthly intervals has been estimated in this form, as well 
as other indices on an experimental basis. 

7- 1 Alternatives on index return modeling 

It is possible to relax the constant returns assumption for indices and estimate 
the corresponding row of the A matrix. 

Another approach is to have a relationship from yields to the expected 
return. A simple one is a spread over the 3 month yield, although several 
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yields can be modeled this way. The coefficients could be estimated by a 
regression. 

Altering the volatility assumptions can be done with the regime switching. 



Bond indices can be modeled in the same manner as equity indices. We 
consider some alternatives to this below. 
8-1 Bond portfolios 

A set of bond prices can be computed and output for each time point in each 
scenario. Prom these total returns can be calculated. Given a set of portfolio 
weights, these returns can be combined into a portfolio return. One or more 
such portfolios can be calculated. 

We can calculate a portfolio return as 



where W{ are the portfolio weights and the Ri are the portfolio returns. 
One assumption on the Wi is a set of constants. Another assumption is 



8.8 Bond Index Modeling 




(8.19) 



8.8. BOND INDEX MODELING 
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that the amount invested in each bond is constant. Interest is then allo- 
cated to that bond or the portfolio. This assumption is more complicated 
to model. When a bond matures it repurchases into the same maturity is 
usually preferable to the portfolio or a fixed new purchase portfolio. 

A modeled portfolio return can be a weighted average of such portfolio 
returns and other equity like indices. 



where now the Ri are some combination of equity and bond portfolio returns. 
One of these can be a "residual" or selected return. The portfolio weights 
against constructed benchmark bond yield portfolios can be estimated by 
regression. 

8-2 Estimation 

This model requires estimation or selection of the bond portfolio weights 
from bonds as well as other portfolio weights. A residual return can be used 
to be estimated from data. 




(8.20) 
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A set of portfolio weights could be selected based on average portfolio 
composition over some period of time. 

8.9 Bond Portfolios Detail 
9-1 Inputs 

The portfolio weights of a set of bond portfolios in terms of primary bonds. 
The maturity dates, spreads and other information on these primary bonds. 

The portfolio weights from each output index in terms of the bond port- 
folios. 

A "residual" return portfolio in terms of its mean, standard deviation of 
return and correlation to other elements of the state vector y. 

The mean could in theory be a function of y or of other variables such as 
yields of bonds (or even portfolios). 

9-2 Process 

The residual return is calculated as part of generating the vector y. 

The bond portfolio returns are generated by first generating prices of the 
bonds, then calculating their returns and then combining these into the bond 
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portfolio returns. 

The residual return and bond portfolio returns are combined using the 
weights for each index. 

The weights to form the bond portfolio returns and the residual return 
might vary with y or other inputs or outputs, but in the simplest case would 
not. 

9-3 Output 

The value of the index. 

8.10 LIBOR modeling 

LIBOR can be modeled as the primary yield curve. Treasuries can then be 
modeled if needed by a linear regression. Alternatively, one can go from such 
a treasury linear regression to LIBOR. 

In practice it is often sufficient to start with the LIBOR-swap curve and 
to modify the ultimate target rate from treasury parameters by increasing it 
modestly, e.g. by 50bp. The generated LIBOR-swap curve can then be used 
to construct a treasury curve by subtracting the initial spread from treasury 
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to LIBOR-swap. 

If LIBOR is used as the benchmark yield curve, LIBOR bond portfolio 
returns will more closely track an investment grade corporate bond portfolio. 
10-1 LIBOR data 

If LIBOR-swap data is provided, the LIBOR-swap rate process can be better 
estimated than simply modifying treasury parameters subjectively 
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Notation 



9.1 Basic parameters and variables 
1-1 n 

The number of individual data series, e.g. S & P 500, Ford, etc. These are 
numbered 1 to n. . 
1-2 T 

The number of time periods, numbered from 1 to T. 
1-3 x 

An n by 1 vector from a multivariate normal distribution. 
1-4 y 

An n by 1 vector of variables transformed from x. These are the observables. 
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1 - 5 S 

Correlation matrix that is n by n. Used to do simulation, and is the theoret- 
ical unobserved true correlation matrix of the x vector. 

1- 6 W 

This -is the estimated correlation matrix. It is used to do the actual simulation 
as if it was the theoretical correlation matrix when the software system is 
used in practice. 

9.2 Returns 

2- 1 Log return 

Given any time series, V u t=l,...,T, such as a price, or other variable which 
is non-negative, we can define its log return as u t 



u t = lnv t — lnv t -\ (9.1) 

for t=2,..,T. We thus "lose" a data point at t=l. 

So for example, y it might be the change in the natural logarithm of the 
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change of Ford's stock price. We calculate it from the observed stock prices, 
P t and P t _i as 

Va = InPt - InPt-i (9.2) 

where i corresponds to Ford, say i=l was Ford. 

To go the other way, to get P t given y it we calculate 

Pt = Pt-ie*" (9.3) 

We have to know P t _ x . If we are simulating we have to know P 0 say and 
then we can simulate for later t using yii,...,j/ir- 
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Chapter 10 

Best modes: RS-VAR 



0-2 Introduction to a VAR with Regime Switching 

The following is given as an introduction to the idea. We shall however, use 
the definitions that appear later to define this. Suppose we have two state 
vectors, y and v. We have the vector of continuous variables 

dv = (b + Av)dt + Gdw (10.1) 

In addition to this is the state index s, which we interpret as a set of a finite 
number of states. We have a mapping from s to b,A and G, assigning for 
each s, b(s), A(s) and G(s). Different values of s can then correspond to 
the different paramters of the univariate models for each index. We have a 
process on s, of the form p(s', s)dt where p(s\ s) is the transition rate per 
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unit time to go from state s to state s\ In addition the transformation 
functions from v to y are indexed by s in addition to their other parameters 
so y = h(v,s,t) is the transformation function from v to y contingent on 
state s. 

0-3 RS-VAR 

Definition 10.1 (Discrete RS-VAR) Suppose we have a discrete variable 
s, that varies over some set, which can be finite. This might be 1 to p or 
0 to p-1. The variable s can be thought of as a state or a regime. The 
transition probability matrix for transition from a state s to s' is given by 
(/>{s\$). (Note that the order can appear differently in some treatments.) 
The standard requirement is that the sum over different final states s' is 1 
for any initial s\ Each element of (j){s\ s) must be non-negative, but zero 
is allowed. There is a mapping from s to n by 1 vectors // 0 (s), and n by n 
matrices A(s) and G(s). The n by 1 vector y follows the process 

y[t) - y[t - 1] = to[s] + K[s]y[t - l})At + E[s}VAt (10.2) 



Ill 

The time interval can vary with time t. It is also possible to allow /lx,K, S to 
depend on t as well as s. 4 

Definition 10.2 (Continuous RS-VAR) The regime s can make a tran- 
sition at any point t in continuous time. The probability (conditional on a 
transition occurring) of a transition from s to s' is ((s\ s). The probability of 
a transition out of state s is r)(s)dt. In this case, we require that C(s, s) = 0. 
We also require that the ((s\ s) are non-negative and sum to one. 
There is a vector of state variables, v, that follows the process 

dv = (b[s] + A[s}v)dt + G[s]dw (10.3) 

where v is an n by 1 vector of variables, b is an n by 1 vector of parameters, 
A is n by n, G is n by k and dw is k by 1. The vector dw is a vector of 
Wiener processes, with mean 0 and variance dt. 4 

Definition 10.3 (Mixed RS-VAR) The state transition can be fixed at 
time nodes and follow a discrete process as in the discrete Regime Switching 
portion of the discrete RS-VAR. From one time node to the next, the process 
is then a continuous-time VAR. 4 
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Definition 10.4 (RS-VAR) A RS-VAR is either a discrete RS-VAR or a 
continuous RS-VAR. 4 

Definition 10.5 (Essential RS-VAR) An RS-VAR that is not a VAR is 
called an Essential RS-VAR or E-RS-VAR. * 
0-4 Discretization Method 

Definition 10.6 (Discretization method) For a continuous-time stochas- 
tic process unless otherwise specified the discretization method shall mean 
the method obtained by restating the process as a discrete time process for 
all variables using the already determined values from the previous time node 
to calculate all new values at the next time node. In terms of equations, all 
variables on the right of the equal sign are understood as from the previous 
time noe and all on the left as the next time node. Note that dt becomes 
Af and Wiener or diffusion type elements, commonly denoted by dw, dz, etc 
become y/At times an appropriate random variable. Those variables may be 
correlated or uncorrelated as the convention, text, or references indicate. In 
general it is to be understood they can be correlated with some correlation 
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matrix. They are to be understood as mean zero, with individual variable 
variances of one and some correlation matrix. That may also be expressed 
in terms of a Cholesky Decomposition of a variance-covariance matrix. In- 
puts and outputs are to be understood as the appropriate set of varibles and 
parameters for the convention used. The discretization method is a fall back 
in case no other method is specified, available, in the prior art or literature. 
Each claim, element of the specification, etc. that can be so understood in- 
cludes this method as within its meaning except where indicated otherwise 
or it is clear by the context. 4 



0-5 FX 

In an FX model, the parameters can vary with regime. This can be done 
with all the parameters in the FX model considered earlier in reference to the 
MFC-FX model. A multiple set of currency exchange rates can be modelled 
with regime switching parameters. 
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0-6 Defining Regimes 

A single discrete index can be used even if we start with multiple discrete 
indices. Suppose there are n indices, and for the i-th index it can have m* 
possible values. These need not be integers but we can always map them 
onto integers. We then have 



possible states. The product is over i=0 to n-1 or i=l to n depending on our 
indexing scheme. Other schemes are possible as well, of course. It may be 
that there are a smaller number of possible states because some combinations 
can't occur. We can still treat those as states but the probability to go to 
them is zero. If we don't start in one, we won't go to one. We could also 
define the probability of exit from such a state to be one. Let M < N be the 
accessible states. We can then index the actual states s from 0 to M-l or 1 
to M or some other set of values, which need not be integers but are more 
convenient that way, especially for computer programming. 



N = Uiirii 



(10.4) 



10.1. RS-VARDMRP 
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10.1 RS-VAR DMRP 



The specific case of the RS-VAR DMRP is given as follows. We have as 
before, the DMRP as: 



We can then make the coefficients functions of a regime index s, where the 
regime index follows a pure regime switching process in discrete or continuous 
time or with discrete time points at which it can change. Here the short rate is 
r = e u . Note one can introduce a time-dependent factor, so that r = e u a(t). 

Individual yields can be modeled through solving the DMRP. One can 
also have processes on "yield residuals" such as an AR(1) process or joint 
VAR on these yield residuals as well. These can be fit to the initial yield 
curve and decay according to such a process, or a straight line decay over 
a finite time interval or some deterministic pattern of decay or that times a 



du = k(6 - u)dt + adz 1 



(10.5) 



d6 = k 2 (02 - u)dt + a 2 dz 2 



(10.6) 
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random residual which follows an ARIMA process or Vector ARIMA for a 
family of yields. 

These methods can be used for any of the RS-VAR's considered here, 
as can the time-dependent multiplier for the short rate, regardless of the 
transformation method or time-dependent coefficients in any of the functions. 
See the prior art for many examples of this. 

4 

We can also imbed this model into a RS-VAR of n-dimensions 

10-2 Auxiliary processing 

After the production of state variable files or data, and/or yield/price curve 
files or data as scenarios or grids or otherwise, additional scenario files and 
other data can be produced. This can use transformations and functions ' 
involving multiple variables at once from all of these, as well as lagged val- 
ues and using Generalized Financial Variables as defined elsewhere in this 
document. 



10.3. PARAMETER ESTIMATION 

10.3 Parameter Estimation 
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This can be done by a variety of means. These include Generalized Method of 
Moments(GMM), Maximum Likelihood (ML), and also the use of judgement, 
stylized facts, and other methods. See the references and the references cited 
therein for the application of these methods in finance and to VAR's, regime 
switching and in some cases RS- VAR's. 

10.4 Other data institution 

An institution's other data can be combined with the data for the stochastic 
process for performing analysis, producing reports, making preparation for 
a transaction, and executing or attempting to execute a transaction. A 
transaction might also include a bid or ask. This may involve transmitting 
data, or performing a transaction or action over or though or by the use of 
the internet or similar network or electronic or digital methods. See also the 
section on individuals, which may include institutional use for its customers, 
prospective customers or individuals it wishes to transact with. These can 
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be combined as well. 

10.5 Other data individual 

The stochastic process can be used by the individual, an adviser, an institu- 
tion for which they are a customer or prospective customer or that desires to 
make some transaction with the individual. This can be combined with data 
for that individual or data that might be used for a representative individual, 
illustration, benchmark or profile of that individual or of one or more cate- 
gories used to prepare the analysis, quote, transaction, etc. for or with the 
individual. This may involve transmitting data, or performing a transaction 
or action over or though or by the use of the internet or similar network or 
electronic or digital methods. 



Chapter 11 

Best Modes: Definitions 



Definition 11.1 (Computer) A device with one or more of the following 
elements 



1. Hard Drive 



2. Memory 



3. Central Processor Unit(s) 



4. Electronic Chip 



5. Pathways for conducting electrons encoded into a physical object. 



6. Cathode Ray Tube 
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7. Printer whose output is controlled by digital media, software, or elec- 
trons 

8. Computer chips such as Intel 80286 architecture based, or its descen- 
dants, or its predecessors, or its competitors. 

9. Computer chips such as Intel 80386, 80486, Pentium, Pentium II, etc. 

10. Devices incorporating these or similar chips. 

11. Devices capable of operating Microsoft operating systems, including 
DOS, Windows 95, 98, 2000, NT, XP, and ones incorporating the code 
or methods or replicating some of the operation of these systems. 

12. Devices capable of operating LINUX, UNIX or similar systems. 

13. Devices capable of operating operating systems used for DEC, HP, 
Apple, IBM or similarly marketed computers from personal computers 
to mainframes. 

14. Devices capable of being operated by similar or competing operating ' 
systems. 
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15. Devices connected to devices operated by the aforesaid or similar op- 
erating systems. 4 

cite some handbooks and product numbers?? Computing device and other 
synonyms shall be understood as referring to this definition. 

Definition 11.2 (Internet) Means for transmitting electronic information 
involving one or more of the following elements 

1. Transmission Control Protocol: TCP 

2. Internet Protocol: IP 

3. Network protocol 

4. Cables to conduct electrons between computers 

5. linked computers 

6. methods to link computers 

7. devices to link computers 
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8. Any network part of which is appropriately documented by the CISCO 
Internetworking Technology Handbook. 

4 

9. Any network capable of performing a task also performed by the fore- 
going. ^ 

See CISCO Documentation Internetworking Technology Handbook, espe- 
cially Chapter 30, Internet Protocols. 

Network and other synonyms shall be understood as referring to this defini- 
tion. 

Definition 11.3 (Equation) In the context of forming computer or net- 
work algorithms, an equation shall have the following meaning. Except where 
an implicit method is indicated an equation shall have the same meaning as 
in a programming language like C, C++, FORTRAN, or higher or lower level 
computer languages. It shall mean the machine or the process as appropri- 
ate to the context of the claims for computer related processes, machines or 
patentable subject matter as in patent office guidelines, 705 type patents, 
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and decisions of the CAFC or other courts. Where an equation is used for an 
implicit algorithm it means the calculation of the disrepancy of the equation 

r 

or deviation as appropriate to an implicit algorithm. In the case of stochas- 
tic differential equations or partial differential equations or other equations 
requiring an algorithm, a standard algorithm or an algorithm in the patent 
shall be understood as appropriate means unless the context clearly indicates 
some other meaning. ^ 

Note 11.1 (Diagrams as Computer or Process Diagrams) In the con- 
text of forming computer algorithms or stating the use of a computer or net- 
work a diagram shall have the following meaning. In such context, a diagram 
or several diagams shall be understood as referring to processes involving 
computers, machines or other patentable processes, machines, methods, etc 
as are common in 705 type patents, patent office guidelines, decisions of the 
CAFC and other courts. ^ 

Definition 11.4 (Function) In the context of forming computer or net- 
work algorithms, a function shall have the following meaning. A function in 
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this context means as in C, C++, etc. a set of operations or a computer 
function or subroutine and the appropriate encoding of it onto a computer. 
In some cases, an implicit algorithm is indicated by the context or the prior 
art of this patent. 4fc 

Note 11.2 (References as Computer Calculations) Where references (i. 
texts like a book or article) are referred to, their use of equations, functions, 
calculations, algorithms, etc. when used as part of the specification or patent 
shall be understood as referring to the definitions of function, equation, cal- 
culation, etc. given here. 4 

Definition 11.5 (Calculation) In the context of forming computer or net- 
work algorithms, calculation means using a machine programmed for that 
purpose including a computer or if appropriate a computer. 4 

Definition 11.6 (Variable Transformations) Given a vector i>,i=0,...,n- 
1, we could make a set of transformations such as 



3. x k = ln(\vi\) 

4. x k = cos(vi) 

5. x k = sin(Vi) 

6. Xk — tan(vi) 

7. x { = Zn(v / H) 

.8. x p = , a real where this is defined. 
9. x p = |^| a , a real where this is defined. 

10. x p = I( Vi > 0) 

11. x p = sgn(vi), where this is the sign of i^, 

12. etc. 

13. recursive functions using polynomial, rational, analytic, indicator, and 
inverse function or relations, integral, differences, weighted sums, etc. 
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Note that in the above, the index in the x vector need not correspond to that 
of the v vector. The repetition of the same index in the above is not relevant. 
It means a separate transformation. The above can be real functions, from 
reals to reals, where defined or from complex to complex where defined. In 
the latter case, there may be poles, branches and a Riemann Surface. We 
can also make transformations involving several of the v's at the same time. 

1. r = v'Qv + 0'v + a 

2. r = e v 'Qv+0'v+a 

3. Multivariate Copula 

4. etc. 

5. recursive functions using polynomial, rational, analytic, indicator, in- 
tegral, differences, weighted sums, etc. 

Here x h etc. r could be an interest rate, default rate, prepayment rate, lapse 
rate, expense, etc. The parameters above can themselves be functions of 
the variables. Complex numbers or transforms to complex variables or other 
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objects such as vectors or matrices of complex numbers, or sequences or series 
of them, etc. are allowed. The variables themselves can be complex even in 
the RS-VAR. 

We can extend the state vector of say a RS-VAR by using lagged elements, 
using an integral or sum of functions of the elements of the RS-VAR, and 
past or current values of those or other variables. A 

Definition 11.7 (GFV: General Financial Variables) AGFVisanyof 
the variables from the following list 

1. Any variable appearing in 

(a) A public report of a company relating to it, such as 

i. Annual Report 

ii. Income Statement 

iii. Statement of Changes in Financial Position 

iv. Notes, Appendix or reference to any such statement. 

v. The same prepared with a different system of accounting in- 
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eluding mark-to-market accounting, mark-to-model account- 
ing, or mixtures of them. 

(b) A filing with 

i. SEC 

ii. Federal Reserve 

iii. State Insurance Commissioner 

iv. IRS 

v. similar entity 

vi. etc. 

vii. The same prepared with a different system of accounting in- 
cluding mark-to-market accounting, mark-to-model account- 
ing, or mixtures of them. 

(c) Any variable prepared in accordance with any of the following 
accounting methods 

i. US GAAP 

ii. Canadian GAAP 
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iii. GAAP of any country or union of countries such as the Eu- 
ropean Union 

iv. International Accounting Standards 

v. US Statutory Accounting based on NAIC,AAA, SOA, CAS, 
CIA, or other actuarial organization or state insurance or 
Canadian federal or provincial Insurance regulation. 

vi. The same for a European country or politicial, economic or 
other division thereof. 

vii. Any of the above using a modification involving mark-to- 
market or mark-to-model or a mixture of them. 

(d) etc. 

2. Managerial Accounting or internal financial analysis. 

3. Financial analysis variables of an external group, such as those created 
in the course of preparing the following. A variable 

(a) Prepared for or by a rating agency 
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(b) Prepared for or by a broker or dealer in securities, derivatives, 
exchange traded instruments, or financial contracts. 

(c) Prepared for or by an accounting or consulting firm. 

(d) Prepared for or by an actuarial firm. 

Specific instances from the list below 

(a) Price 

(b) shadow price 

(c) market price 

(d) model price 

(e) theoretical price 

(f) price curve varying a parameter or other prices 

(g) cash flow 

(h) cash payment 

(i) valuation 
(j) reserve 
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(k) capital 
(1) percentile 
(m) conditional tail expectation 

(n) modified conditional tail expectation (as in US American Academy 
of Actuaries C3 Phase II, i.e. cut off at zero in gain for any sce- 
nario) 

(o) required capital 

(p) required reserve 

(q) policy holder behavior 

(r) withdrawal rate 

(s) new business rates 

(t) policy holder elections 

(u) contract behavior variables 

(v) business behavior variables 

(w) corporate behavior variables 
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Additional specific instances (counter too large for one list) 

(a) market behavior variables 

(b) government behavior variables 

(c) tax rate(s) 

(d) deductions and exclusions 

(e) tax payer utilization of deductions and exclusions 

(f) amount due 

(g) settlement amount 

(h) mortality 

(i) lapse 

(j) prepayment 
(k) exchange rate or rates 

(1) default rate or rates 
(m) reserve against default 
(n) credit migration 
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(o) other migration variable (prepay, lapse, expense etc.) 
(p) interest 
(q) accrual 

(r) a value used in computing tax or from a computation of tax 
whether a tax liability, accrual, cash amount due, penalty, etc. 

Mortality, lapse, prepayment, etc. might be as rates or quantities. 
Multiple groups or types of any of these or other variables are included 
in this definition. 

In a legal context, any variable relating to or used in the determination 
of any of the following: 

(a) Valuation of any amount in a legal framework 

(b) In Contracts, settlement of contracts, damages or restitution 

(c) In Torts.. 

(d) Any quantity, numeric or otherwise subject to mathematical anal- 
ysis or projection that is part of the law of remedies, damages, 
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restitution, etc. 

(e) Real property, mortgages, debt, liens, etc. 

(f) Quantity in a statute 4 

6. Any intermeidate variable or value used or useful to calculate any of 
the above. 

7. Any input variable to any calculation of any of the above variables. 

8. Any output variable from any calculation used to calculate or ulti- 
mately calculate any of the above variables. 

Definition 11.8 (Financial Variable) We shall consider the financial vari- 
able and GFV to be the same. ^ 

Definition 11.9 (Exponential GFV) An exponential GFV is one where 
to calculate at least one of the GFV's an exponential function or a numer- 
ical algorithm equivalent to such a function call in a standard language is 
performed. A 
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Definition 11.10 (Short Term Interest Rate Exponential GFV) A short- 
term interest rate exponential GFV is one where to calculate the short term 
interest rate, a call to an expnential function is made, or a numerical algo- 
rithm equivalent to such a call in a standard language is performed. A 

Note 11.3 (Algorithm Version) Where an algorithm is not specified for 
an equation or diagram the following are to be understood as one of the 
algorithms for that type of calculation. Examples are 

1. Discretization 

2. Using a "closed form formula" (CFF) 

3. Using a CFF as an approximation 

4. The corresponding integral/sum form and its corresponding discretiza- 
tion. 

5. Taking the corresponding integral/sum form and using approxima- 
tions to "Green's Functions" or Fundamental Solutions with or without 
boundary conditions or special side conditions. 
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6. Using numerical algorithms such as in SIAM journals and publications 
or other standard references. 4 

The algorithms in turn are to be understood as referring to patentable subject 
matter, i.e. process or machine or method such as on a digital computer 
programmed with these algorithms. 

Definition 11.11 (Time Loop Initiation Step) Information is reset to , 
the initial values as appropriate. This includes resetting counters, memory, 
arrays, registers, etc. 

Definition 11.12 (Time Loop Recursion Step) Information from the prior 
time node is updated to the next time node. 

Definition 11.13 (Scenario Initiation Step) Information 



Definition 11.14 (Next Scenario Step) Information 



137 

Definition 11.15 (Scenario Combinations) Information from different sce- 
narios is combined in some applications. For example, a price in some calcu- 
lations is the average of the discounted cash flows in different scenarios where 
those discounted cash flows are calculated using the appropriate variables for 
that scenario, such as cash flows, discount rates, default rates, lapse and pre- 
payment rates or the like, expenses, taxes, commission, fees, penalties, etc. 

Definition 11.16 (Recursive Monte-Carlo) A Monte Carlo that calls a 
monte carlo. Also called a Monte Carlo within a Monte Carlo. Either or both 
can use QRMC/LDS or convention random number generation methods. 
This could be called a simulation within a simulation. Multiple such nested 
calls are allowed. The word simulation and its synonyms shall allow for this 
recursive structure in this document. A 

Definition 11.17 (GFV Simulation) A GFV Simulation Process: A sim- 
ulation process to simulate a GFV. It can involve input of random variables 
to calculate the GFV or their probabilities or other measures. This is a set 
of steps to follow. ▲ 
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Definition 11.18 (GFV-SS: GFV Stochastic Simulator) A GFV Stochas- 
tic Simulator (GFV-SS) is defined as: 

A GFV Simulation Process based on a stochastic process such as a VAR 
or RS-VAR. 

Definition 11.19 (RS-VAR GFV-SS: RS-VAR GFV Stochastic Simulator) 

A RS-VAR GFV Stochastic Simulator (GFV-SS) is defined as: a GFV 
Stochastic simulator that takes output from or uses a RS-VAR. The RS- 
VAR can be run prior to running another module of the GFV-SS or during 
it. The RS-VAR can produce a file or data in memory. It can be run be- 
fore the remaining modules or another module. In a subsequent module, 
the RS-VAR output can be read from a file or calculated during the other 
module. Or the other module could use a formula based on the RS-VAR, or 
calculating a probability or other measure. A characteristic function, Green's 
function, state price, martingale, or other function can be a module or part 
of a module. The RS-VAR or other modules can use the so-called P (ob- 
jective) or Q (risk-neutral) measures or other measures. It can use multiple 
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measures, even inconsistent ones in the same or different modules. 4 

Definition 11.20 (RS-VAR GFV Stochastic Simulator Machine) A com- 
puter programmed with an RS-VAR GFV Stochastic Simulator. * 

Definition 11.21 (RS-VAR GFV Stochastic Simulator Article of Manufacture) 

A physically tangible item manufactured or altered by the use of a RS-VAR 
GFV Stochastic Simulator Process or a RS-VAR GFV Stochastic Simulator 
Machine. A 

Definition 11.22 (Financial Product) 1. Bond 

2. Mortgage 

3. CMO 

4. Bank account 

5. Exchange traded instrument. 

6. A financial contract. 

7. An ISDA contract. 
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8. A regulated financial contract. 

9. An insurance contract. 

10. Annuity. 

11. Universal Life. 

12. Whole Life. 

13. Guarantee. 

14. Rider on a policy or contract. 

15. etc. * 

Definition 11.23 (Transaction End-Use (TEU)) An action consisting 
of one or more of the following actions or the actualization of a result in this 
list. 

1. Preparing data for use in preparing any report listed in the GFV defi- ' 
nition. 
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2. The value of a GFV in a report or file. 

3. A transaction involving the purchase or sale of a financial product. 

4. A payment related to a financial product. 

5. A bid or ask quotation of a financial product. 

6. A transaction of a sale, lease, loan, repurchase agreement, formation 
or sale of a collateralized obligation, etc. 

7. A quotation for the same. 

8. Preparing one or both of a bid or ask. 

9. Transmitting the bid or ask. 

10. Accepting a bid or ask. 

11. Preparing an indication or appraisal. 

12. A transaction subject to an inequality restriction such as a purchase at 
a lower value than a calculated quantity 
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13. A transaction subject to a measure, such as a purchase based on a 
measure of the deviation above a calculated quantity. 

14. Valuing inventory of financial instruments or valuing financial con- 
tracts. 

15. Valuing liabilities or assets. 

16. Publication of the values or intermediate values from or associated with 
the above. 

17. Making a purchase, offer, sale, or bid on the internet. 

18. Making a purchase, offer, sale, or bid by electronic means. 

19. Preparing a report as part of a service of consultation or advisory work. 

20. Preparing a report as part of a financial service. 

21. Preparing a report as part of a regulated financial service. 

22. Advising a purchase, sale, or exercise of an option as part of a financial 
service. 
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23. Determing a dividend payment based on a financial report. 

24. The paying of a dividend determined in said manner. 

25. Determining an amount of a payment based on a financial report. 

26. The paying of said amount. 

27. Determing a dividend payment based on a financial report that involves 
quantities reflecting a simulation of future financial condition. 

28. Determing a dividend payment based on a financial report that involves 
quantities reflecting a simulation based on an Essential VAR. 

29. The making of said payment. 

30. Determing any payment based on a financial report that involves quan- 
tities reflecting a simulation based on an Essential VAR. 

31. The making of said payment. 

32. Determining a credit rating or regulatory condition based on a simula- 
tion of an Essential RS-VAR. 
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33. A purchase or sale of bonds, currency or other securities by a central 
bank. 

34. A purchase or sale of bonds, currency or other securities by a central - 
bank to change a money supply. 

35. A purchase or sale of bonds, currency or other securities by a central 
bank resulting in a change in money supply. 

36. Altering an interest rate or credited rate or terms of interest. 

37. The preceding done by a central bank. 

38. Providing a report or rating to a counter-party or customer or a po- 
tential counter-party or customer. 

39. Providing a report to a government regulator, exchange or private reg- 
ulatory body. 

40. Preparing such a report. 

41. Running a data processing system for the above. 
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42. Running a "data Processing System for Hub and Spoke Financial Ser- 
vices Configuration" for the above. 



43. Attempting any of the above. 



44. Preparation for the above. 



45. Using an Essential RS-VAR in any of the above. 



46. Using a GFV calculated from an Essential RS-VAR in any of the above. 



47. Using a stochastic simulator in any of the above. 



48. Any of the above using or that used a stochastic simulator. 



49. Any of the above using or that used an Essential RS-VAR, or a GFV 
computed by using an Essential RS-VAR. 



50. Any of the above transactions using a computer. 



51. Any of the above transactions using data transmission on internet. 



52. Executing a transaction such as the above on the internet. 
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53. Executing a transaction such as the above on a network. 

Definition 11.24 (Transaction End Use Quantity (TEUQ)) A trans- 
action end use quantity is any quantity appearing in the transaction end use. * 
This includes any variable from the following list 

1. A price appearing in a TEU. 

2. A quantity appearing in a TEU. 

3. A variable relating to a TEU. 

4. A variable relating to a financial product in a TEU. A 

A computer or the internet can be used to compute a TEU quantity 
and that TEU quantity encoded onto a physical medium. Such a physical 
medium might be a computer readable medium or memory. 

Elements of the above may be transacted in different places. For example 
a simulation of future financial condition might be made in one country based 
on financial contracts in another country and a payment made in a third 
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country. This might be reported in another country. Dividend payments 
might then be calculated in another country, reported in yet another country 
and finally paid in another country. Rating agency reports might similarly 
be prepared somewhere else and reported somewhere else. These might then 
be used by a counter-party in a derivatives contract, a customer, an advisor, 
or a financial service firm for advising its customers or managing products 
or portfolios for them. 

Definition 11.25 (End Use Entity (EUE)) An end use entity is an or- 
ganization or legal entity that uses TEU's, TEUQ's, or the methods, pro- 
cesses, articles of manufacture, machines of this patent. In addition it may 
be an object that is created, modified, maintained, altered, reported on or 
used by such an EUE or a counterparty or customer or payment recipient of 
an EUE. 

1. Database 



2. Database of Financial Products 
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3. Database of Financial Product data 

4. Database of TEUQ's. 

5. Database containing TEUQ's. 

6. Database prepared with an Essential RS-VAR. 

7. Financial Services Firm. 

8. Line of Business. 

9. Law firm. 

10. Securities class action litigation valuation data. 

11. Class a 

12. Rating Agency. 

13. Ratings report, database or software program accessing it. 

14. A portfolio of financial products. 

15. A trading desk. 



16. A trading desk's positions. 

17. A trader's positions. 

18. Central Bank. 

19. Portfolio. 

20. Pension Fund. 

21. Pension Portfolio. 

22. Money Supply. 

23. Portfolio of Debt. 

24. Bank. 

25. Insurance Company. 

26. Derivatives Portfolio. 

27. ISDA contract counterparty. 

28. Portfolio containing ISDA contracts. 



150 CHAPTER 11. BEST MODES: DEFINITIONS 

29. Portfolio containing Exchange Traded Contracts. 

30. A physical object encoded with any of the above that can be encoded 
onto a physical object. 

31. A machine readable version of said physical object. 

32. "Hub and spokes financial services configuration." 4 
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